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Abstract 

Background: The family of TGF-/J ligands is large and its members are involved in many different signaling 
processes. These signaling processes strongly differ in type with TGF-/3 ligands eliciting both sustained or transient 
responses. Members of the TGF-/3 family can also act as morphogen and cellular responses would then be 
expected to provide a direct read-out of the extracellular ligand concentration. A number of different models have 
been proposed to reconcile these different behaviours. We were interested to define the set of minimal 
modifications that are required to change the type of signal processing in the TGF-/3 signaling network. 

Results: To define the key aspects for signaling plasticity we focused on the core of the TGF-/J signaling network. 
With the help of a parameter screen we identified ranges of kinetic parameters and protein concentrations that 
give rise to transient, sustained, or oscillatory responses to constant stimuli, as well as those parameter ranges that 
enable a proportional response to time-varying ligand concentrations (as expected in the read-out of 
morphogens). A combination of a strong negative feedback and fast shuttling to the nucleus biases signaling to a 
transient rather than a sustained response, while oscillations were obtained if ligand binding to the receptor is 
weak and the turn-over of the l-Smad is fast. A proportional read-out required inefficient receptor activation in 
addition to a low affinity of receptor-ligand binding. We find that targeted modification of single parameters 
suffices to alter the response type. The intensity of a constant signal (i.e. the ligand concentration), on the other 
hand, affected only the strength but not the type of the response. 

Conclusions: The architecture of the TGF-/3 pathway enables the observed signaling plasticity. The observed range 
of signaling outputs to TGF-/J ligand in different cell types and under different conditions can be explained with 
differences in cellular protein concentrations and with changes in effective rate constants due to cross-talk with 
other signaling pathways. It will be interesting to uncover the exact cellular differences as well as the details of the 
cross-talks in future work. 
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Background 

Transforming growth factor beta (TGF-/3) signaling has 
been implicated as an important regulator of almost all 
major cell behaviors, including proliferation, differentia- 
tion, cell death, and motility [1]. Which response is 
induced or repressed depends on the cell type and con- 
text in which the signal is received. 

The complexity of the biological outcomes elicited by 
TGF-/3 stands in stark contrast to the apparent simpli- 
city of the signaling cascade. In response to TGF-/3, type 
1 (ALKs 1-7 in humans) and type 2 receptors (ActR- 
IIA, ActR-IIB, BMPR-II, AMHR-II and TbR-II, in 
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humans) form complexes and the constitutively active 
type 2 serine/threonine kinase phosphorylates the type 1 
receptor. The activated type 1 receptor transduces the 
signal into the cell by phosphorylating the regulatory 
Smads (R-Smad: Smad 2 and 3 in case of the TGF-/3 
subfamily, and Smad 1, 5 and 8 for the BMP subfamily). 
Once activated R-Smads form homomeric complexes 
and heteromeric complexes with the common Smad, 
Co-Smad (Smad 4) [2]. Smads continuously shuttle 
between nucleus and cytoplasm [3]. TGF-/3 signaling 
biases Smad localisation to the nucleus [4] where Smad 
complexes associate with chromatin and regulate the 
transcription of hundreds of genes [5]. Signal termina- 
tion is achieved through continuous dephosphorylation 
of the R-Smad (mainly in the nucleus [3]) and induction 
of inhibitory Smads (I-Smad: Smad 6 for the BMP sub- 
family, and Smad7 for the TGF-/3 subfamily). I-Smads 
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act through diverse mechanisms: by targeting active 
receptor for proteasomal degradation [6,7], inducing 
receptor dephosphorylation [8] and competing with R- 
Smad for the receptor binding site [9]. Rapid shuttling 
and inactivation enables a continuous sensing of the 
extracellular ligand concentrations [3]. This is likely to 
be particular important when members of the TGF-/3 
ligand family acts as morphogen and determine cell-fate 
in a concentration-dependent manner. 

Beyond the core components of this signaling pathway 
many other factors modulate the signal and thereby 
contribute to the versality of the response. At the mem- 
brane level, the access to receptor is controlled by solu- 
ble proteins that sequester TGF-/3 ligand (i.e. decorin) 
[10], and by membrane-bound co-receptors that pro- 
mote binding (i.e. betaglycan) [11]. The receptor activity 
is further regulated by several receptor internalization 
routes [12], and by receptor turnover. Intracellularly, 
many processes require auxiliary proteins (i.e. SARA for 
the binding of R-Smad to the receptor and Schnurri for 
the binding of the R-Smad/Co-Smad complex to the 
DNA binding element) [2,13]. The restriction of those 
auxiliary factors to specific cell-types will make the 
response cell context dependent [14]. Diversity can also 
be generated by the huge number of different possible 
combinations of type 1 and type 2 receptors [2] and the 
multiple crosstalks of the TGF-/3 signaling cascade with 
other pathways. One example of regulation by cross-talk 
is the phosphorylation of R-Smads in the linker region 
by Ras-activated MAPK [15], calcium-calmodulin- 
dependent protein kinase II [16] or CDKs [17]. Phos- 
phorylation reduces the transcriptional activity of the R- 
Smad [18]. 

Several mathematical models have been developed to 
gain further insights into the complex TGF-/3-dependent 
signaling network [19]. An early model by Clarke and 
co-workers (2006) [20] focused on the nuclear accumu- 
lation of Smad complexes. Their conclusion on the cen- 
tral role of the imbalance between R-Smad 
phosphorylation and dephosphorylation rates were con- 
firmed by a more detailed model by Schmierer et al. 
(2008) [3]. Experiments suggest that the duration of the 
response to a ligand stimulation strongly impacts on the 
cellular response. Thus epithelial cells that elicit sus- 
tained nuclear Smad complex accumulation respond to 
TGF-/3 with cell growth arrest, whereas pancreatic 
tumor cells that elicit a transient response continue pro- 
liferating (while keeping other TGF-/3 induced beha- 
viors) [21]. Much theoretical work therefore focused on 
how sustained, transient, or switch-like responses could 
be obtained by adjusting the receptor dynamics, ligand 
depletion, and the I-Smad dependent negative feedback. 
Melke et al. (2006) [22] focused on the potential role of 
I-Smads in generating transient responses while Vilar et 



al. (2006) focused on the receptor dynamics to explain 
the occurrence of both transient and sustained 
responses. Zi et al. (2007) [23] included a simple model 
of the Smad dynamics and highlighted the importance 
of the balance between clathrin-dependent endocytosis 
and non-clathrin mediated endocytosis. All pathway ele- 
ments were finally brought together by Chung et al. 
(2009) [24] in a more comprehensive model, used to 
examine the contradictory roles of TGF-/3 in cancer pro- 
gression. Lately Zi et al. (2011) [25] published a study 
that highlights the potential of TGF-/3 ligand depletion 
in converting short-term graded signaling responses into 
long-term switch-like responses. Unlike for other path- 
ways oscillations have not yet reported for the TGF-/3 
signaling pathway [26,27]. TGF-/3 type ligands are also 
acting as morphogens, and the response to these 
appears to be proportional. Recently, Paulsen and co- 
workers published a study on the impact of synexpres- 
sion of the feedback inhibitors BAMBI, Smad6, and 
Smad7 on the read-out of morphogen gradients during 
embryogenesis [28]. 

While the many published studies explain the different 
behaviours for the different situations for which they are 
observed and highlight the many mechanisms that 
enable the different response types it remains largely 
unclear how easily the response type can be changed. 
We wondered how the TGF-/3 signaling pathway accom- 
plishes the flexibility in its responses and which and 
how many parameters have to be altered for cells to 
respond differently. To efficiently explore the canonical 
response we focused on the core signaling architecture, 
and did not consider the detailed receptor dynamics and 
cross-talks in the model; they are included indirectly 
through the parameters that they modulate. We 
explored the response types and in particular changes in 
the response type as we explored the parameter values 
within biologically meaningful ranges. We find that rela- 
tively small changes in single parameters can alter the 
response. Cellular protein concentrations are a particu- 
lar powerful point of control and this explains how dif- 
ferent cell types can show different responses. 
Importantly we also identify key parameters that affect 
the response and we can relate these to observed points 
of cross-talk between signaling pathways. The particular 
architecture of the TGF-/3 network thus allows for the 
great flexibility in the response. 

Methods 

The model 

Several models for the TGF-/3 signaling network have 
been developed that focus on different aspects of the 
TGF-/3 signaling network, i.e. the receptor dynamics 
[23,29], the shuttling between the cytoplasm and the 
nucleus [3], and the negative feedback via the I-Smad 



Celliere et al. BMC Systems Biology 201 1, 5:184 
http://www.biomedcentral.eom/1 752-0509/5/1 84 



Page 3 of 1 3 



(Smad7/Dad) [22]. These different aspects have lately 
been combined in a model that addresses differences in 
TGF-/3 signaling between normal and cancerous cells 
[24]. The models of the TGF-/3 signaling pathway 
showed that stimulation could result in either transient 
and sustained responses dependent on the choice of 
parameters [3,22-25,29]. Transient responses could be 
obtained through complex receptor dynamic [29], the I- 
Smad-mediated negative feedback [3,22], or ligand 
depletion [25]. Negative feedbacks can in principle also 
give rise to oscillatory behaviour. We wondered whether 
all three qualitative behaviours (sustained, transient, or 
oscillatory response) could be obtained already with the 
most simple intracellular feedback mechanism, and how 
these behaviours would depend on the parameters. 
Since the more complex interactions (that we ignore) 
effectively modulate the parameter values in our model 
an in-depth understanding of the parameter dependen- 
cies in the simple model should also enable a better 
understanding of the complex network interactions that 
are found in the cell. The different response types can 
also (trivially) be obtained by modulating the protein 
concentrations accordingly. We, however, keep the con- 
centrations of receptors, ligand, R-Smad and Co-Smad 
constant and thus include these effects only indirectly as 
changes in the effective binding rates. 

Accordingly, we formulated a detailed model of TGF- 
P signaling that focused on the negative feedback, but 
did not include any complex receptor dynamics as these 
require changes in the receptor and ligand concentra- 
tions. Our model describes the dynamics of TGF-/3 
ligand (TGF-/3), receptor (TGFfiR), regulatory R-Smads 
(denoted simply Smad), Co-Smads, I-Smads, their com- 
plexes as well as the expression intermediates of the I- 
Smad. Importantly, we include two compartments, the 
nucleus and the cytoplasm, and the Smad and Co-Smad 
complexes can shuttle between the two compartments 
as first described in [3]. The regulatory interactions are 
summarized in Figure 1 (a SBML file is provided in 
Additional file 1). Thus the ligand TGF-/3 reversibly 
binds to the TGF-/3 receptor (reactions 1 and 2 in Fig- 
ure 1), which is then phosphorylated to become fully 
active (3 and 4). The active receptor induces phosphory- 
lation of R-Smad (7), which in turn can reversibly 
dimerize or form a complex with Co-Smad (10 and 11). 
Those two reactions can take place either in the cyto- 
plasm or in the nucleus and the five species Smad, 
phosphorylated Smad, Co-Smad, homodimers and het- 
erodimers can shuttle from the cytoplasm to the nucleus 
and back (8, 9 and 12). Nuclear Smad/Co-Smadf com- 
plexes act as transcription factors and trigger the tran- 
scription of I-Smad mRNA in the nucleus (14 and 15). 
The I-Smad mRNA then shuttles to the cytoplasm (16), 
where it can be degraded (17) or translated into I-Smad 



(18) . I-Smad mediates a negative feedback by sequester- 
ing the active receptor (5 and 6) and can be degraded 

(19) . The response to a stimulus by TGF-/3 ligand is a 
change in the transcriptional activity, monitored as the 
nuclear concentration of Smad/Co-Smad complexes. 

We translated those interactions into sets of ODEs 
using the law of mass action where appropriate. To 
reduce the complexity of the model we also employed 
Hill functions to describe the regulation by cooperative 
interactions. To efficiently investigate the impact of 
changes in total concentration of receptors, R-Smad, 
and Co-Smad we used a total concentration rather than 
production and degradation rates for these species. 

To respond to TGF-/3 cells must be able to detect 
changes in the ligand concentration and convert the dif- 
ferences into different transcriptional responses. Tran- 
scriptional activity is determined by the concentration of 
transcription factors in the nucleus. We therefore moni- 
tor the nuclear concentration of R-Smad/Co-Smad com- 
plexes as a measure of transcriptional activity, in 
response to a change in the extracellular TGF-/3 
concentration. 

Parameter screening and simulations 

We are interested in the signaling capacity of the TGF-/3 
pathway within its physiological limits. These physiologi- 
cal limits are set by the plausible range that the para- 
meter values can take. We established a likely range for 
each parameter value based on available data and esti- 
mates (Additional file 2, Table SI and Table S2). While 
previous measurements and estimates are necessarily of 
limited accuracy and differences are likely to exist 
between different cells and different cell types 
[3,22-24,29] we expect that basing ourselves on the 
available data will not too much distort the ranges that 
we screen. Most parameters were varied over 3 or 4 
orders of magnitude, centered around the mean of 
values found in the literature. Since there are no good 
estimates for the I-Smad expression rates kl4 and kl5 
were varied over 5 orders of magnitude. The rates of 
phosphorylation and dephosphorylation (k7 and kl3) 
were varied only over two orders of magnitude because 
a large fraction of the simulations failed when these rate 
constants were varied over a wider range. To avoid a 
bias to the few parameter sets that do not lead to 
extreme dynamics we had to constrain these two para- 
meters to only vary over two orders of magnitude. To 
determine the possible range of pathway responses to a 
defined stimulus, we carried out 10 6 independent simu- 
lations with parameter values randomly picked from a 
uniform logarithmic distribution of parameter values 
within the set ranges (as discussed in Geier et al. [30]) 
and compared the predicted nuclear concentration of R- 
Smad/Co-Smad complexes in response to the ligand 
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Nucleus 



Figure 1 A simple model of TGF-/3 signaling with l-Smad mediated negative feedback. The ligand TGF-/3 reversibly binds to the TGF-/3 
receptor (reactions 1 and 2), which is then phosphorylated to become fully active (3 and 4). The active receptor induces phosphorylation of R- 
Smad (denoted simply Smad)(7), which in turn can reversibly dimerize or form a complex with Co-Smad (10 and 1 1). Those two reactions can 
take place either in the cytoplasm or in the nucleus and the five species Smad, phosphorylated Smad, Co-Smad, homodimers and heterodimers 
can shuttle from the cytoplasm to the nucleus and back (8, 9 and 12). Nuclear Smad/Co-Smadf complexes act as transcription factors and trigger 
the transcription of l-Smad mRNA in the nucleus (14 and 15). The l-Smad mRNA then shuttles to the cytoplasm (16), where it can be degraded 
(17) or translated into l-Smad (18). l-Smad mediates a negative feedback by sequestering the active receptor (5 and 6) and can be degraded 
(19). 



stimulus. In a first step, we let the system equilibrate for 
1 hour with almost no ligand (concentration of 10" 6 pM 
to avoid failure of the solver) and initial cellular concen- 
trations TGF/3R = 1 nM, Smad = 60 nM and Co-Smad 
= 100 nM. We then used the steady-state value of the 
first step and solved the simulations for 10 hours with a 
constant ligand concentration of 200 pM. Using 
MATLAB's odel5 s routine the 10 6 simulations took in 
total approximately 140 hours of CPU time. 

Criteria to define the different TGF-/J signaling responses 

In response to ligand exposure we observed five different 
qualitative responses, i.e. unresponsive, sustained, transi- 
ent, dampened oscillatory or sustained oscillatory 
responses (Figure 2). Additional file 3, Fig. SI, S2, S3, and 



S4 show the evolution of the concentration of each species 
over time in a representative transient and a representative 
sustained response (Additional file 2, Table S3). To define 
the parameter dependency of the different response types 
we made the following definitions: We speak of unrespon- 
siveness if the concentration of nuclear R-Smad/Co-Smad 
complexes remains below a chosen threshold 6 within ten 
hours of stimulation. Accordingly we speak of responsive- 
ness if the concentration exceeds the threshold concentra- 
tion ft and here we distinguished four distinct behaviours, 
inspired by the work of Ma et al.[31] and based on the 
subsequent dynamics: 
7. Sustained response 

After the initial peak the response must retain at least 
90% of its maximal value (called Opeak). To exclude 
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Figure 2 Criteria to define the different TGF-/3 signaling responses. All responses must exceed a threshold concentration of 10 pM initially 
to be considered responsive. (A) Sustained responses: The response must reach 90% of the maximal value Opeak within 2 hours and retain 90% 
of this maximal value until the end of 10 hours simulation. (B) Transient responses: The response must exceed 10 pM and subsequently return 
to less than 10% of the highest value Opeak within 2 hours of stimulation. The final value Oend must be lower than 1 pM. (C-D) Oscillations: 
After the initial >10 pM peak at least four further peaks must exceed 1 pM in amplitude. Depending on whether the fifth amplitude is less or 
higher than half the second amplitude we distinguish (C) sustained and and (D) dampened oscillations respectively. 



slowly increasing responses we further require that 90% 
of the peak value Opeak is reached within less than 
7200 s (2 hours). 

2. Transient response 

After the initial peak the response must drop to levels 
lower than 10% of its peak value Opeak within less than 
two hours and the final value (after 10 hours) (called 
Oend) must be lower than 0.1 x 9. 

3. Oscillations 

After the initial peak the amplitude (difference between 
the local maximum and local minimum) must exceed 
0.1 x 0 at least 4 times. 

3.1 Dampened oscillations The fifth amplitude must be 
less than half the second amplitude. 

3.2 Sustained oscillations The fifth amplitude must be 
higher than half the second amplitude. 

We characterized the long-term behaviour of oscilla- 
tions based on the relative amplitudes of the second and 
fifth peak because the first peak can be particularly high 
(Figure 2D), and most dampened simulations have no 
more than five peaks. 

Quantitative data on the physiological concentrations 
of the cellular proteins and the transcription factor com- 
plex (nuclear Smad/CoSmad complex) do not exist, and 



we therefore had to set our detection threshold arbitra- 
rily to 6 = 10 pM when analysing a unique constant sti- 
mulus with 200 pM TGF-/3 ligand. When the response 
to several ligand concentrations or with several protein 
concentrations was studied we used the maximal 
response value as 9. Simulations were run for 10 hours. 
In case if oscillations, if the amplitude of oscillations 
was still larger than 0.1 x 6 after 10 hours, the simula- 
tion was continued until the oscillations vanished, but 
for a maximum 100 hours. In this way we avoid any 
impact of period length on the classification of oscilla- 
tions, and the length of the period indeed does not bias 
our characterisation of oscillations to dampened or sus- 
tained oscillatory behaviour (Additional file 3, Fig. S5C). 
The time thresholds 2 hours and 10 hours were chosen 
based on experimental data [21]. 

Results and Discussion 

In response to a sustained stimulus (200 pM ligand) our 
simple model for TGF-/3 signaling can give rise to sus- 
tained (Figure 2A), transient (Figure 2B), or oscillatory 
(Figure 2C, D) responses. The sustained/transient dis- 
tinction is particularly relevant, as it has been shown 
that those two qualitative behaviors are related to the 
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growth inhibitory effect of TGF-/3 [21]. To better under- 
stand the conditions for these different behaviours we 
sought to identify parameter families that would give 
rise to a certain response type. We hoped that a com- 
parison of those families would reveal the critical para- 
meters that determine the response type. To that end 
we screened a large number of parameter sets and clas- 
sified them according to their responses as described in 
detail in the Materials and Methods section. 

Parameter-dependent distinct qualitative responses 

Our criteria in Figure 2 are very strict (i.e. the speed of 
responses is an important criterium) and there is a wide 
undefined range between sustained and transient 
responses. As a consequence most parameter sets (57.5%) 
do not fall into any of the defined categories (Additional 
file 3, Fig. S5A, black sets). Of those that can be classified 
most (25.5% of the parameter sets tested) led to no 
response (Additional file 3, Fig. S5A, marked in grey). 
Among the "responsive" parameter sets most lead to sus- 
tained responses (14.8% of the parameter sets tested) 
(Additional file 3, Fig. S5A, marked in red) while transient 
responses are observed less frequently (2.2%of the para- 
meter sets tested) (Additional file 3, Fig. S5A, marked in 
green). All three behaviors have previously been observed 
in various models of TGF-/3 signaling. We find that in 
addition in a minority of cases (306 simulations, i.e 
0.046%) also oscillatory responses can be produced (Addi- 
tional file 3, Fig. S5B). Even though the number of sets 
that give rise to oscillations in the concentration of nuclear 
transcription factor complexes is small, these may occupy 
a sufficiently dense subspace in the parameter space to be 
physiologically relevant. The oscillations can either be sus- 
tained or dampened, depending on how fast their ampli- 
tude decays (Figure 2C, D). As expected sustained 
oscillations have a larger number of peaks (Additional file 
3, Fig. S5B). While the period of the oscillations is not 
biased to dampened or sustained oscillatory behaviour 
(Additional file 3, Fig. S5C), the duration (the time until 
oscillations vanish) depends on both the number of peaks 
and the duration that tends to be higher for sustained 
oscillations (Additional file 3, Fig. S5D). 

Oscillatory behavior has been reported for a number 
of other signaling pathways (i.e. the ERK cascade [32]), 
but so far no experimental evidence exists for oscilla- 
tions in the TGF-/3 pathway. However, standard bio- 
chemical experiments average over a large number of 
non-synchronized cells. If the nuclear concentration of 
transcription factor indeed oscillated, only sophisticated 
single-cell assays would reveal these. 

The impact of kinetic parameters on the response type 

We wondered which kinetic parameters would be criti- 
cal for the different response types. Our sampling space 



is huge (23 parameters, with most of them sampled over 
4 orders of magnitude, Additional file 2, Table SI) and 
we looked for parameters that would be constrained in 
the different response types. In Figure 3 we plot the 
sampled ranges in grey, and the parameter ranges that 
correspond to the different response types in colours. 
Since we are sampling from a uniform logarithmic dis- 
tribution parameters that are not affecting the response 
type should remain uniformly logarithmically distributed 
in the parameter subsets. In Figure 3A we compare the 
parameter ranges of sustained (red) and transient 
(green) responses. We notice that whereas some para- 
meter values remain (almost) uniformly distributed, 
others are constrained. Constrained parameters include 
the rates that describe the I-Smad dependent negative 
feedback loop (parameters £5, £6, £14, £15, £17, £18, 
and £19), the shuttling rate between cytoplasm and 
nucleus (£8), the dynamics of the Smad homo- and het- 
erodimer formation/dissolution(£lO, and £11), and the 
dephosphorylation of Smad (£13). Figure 4 A shows the 
clear segregation of the "sustained" (red) and "transient" 
(transient) parameter sets in a plane spanned by the 
parameters that determine the strength of the negative 
feedback ((£14 x £18)/(£15 x £17 x £19) x k5/k6) and 
the speed of Smad dephosphorylation, £8 x (/ell x £13)/ 
£10. To favour transient responses over the sustained 
responses the I-Smad dependent negative feedback must 
be strong and dephosphorylation of Smad must be fast. 
The need for rapid dephosphorylation likely arises also 
because of our requirement that adaptation must hap- 
pen within 2 hours. We notice that the size of the para- 
meter set that permits transient responses is 
considerably smaller than the parameter set that permits 
sustained responses. However, transient responses can 
also result from degradation of core signaling compo- 
nents and ligand which is not considered here. 

Transient and oscillatory responses are similar in that 
the response must decay quickly in spite of the continu- 
ous presence of ligand. A similar comparison of the 
parameter ranges that permit transient (green) or oscil- 
latory (blue and magenta) responses (Figure 3B) indeed 
reveals that similar restrictions apply (i.e. large shuttling 
rate between cytoplasm and nucleus, £8, and strong 
negative feedback, £14, £15, £17, £18, and £19). How- 
ever, in case of oscillations the response restarts and in 
addition we indeed notice a strong restriction of the 
rate of ligand-receptor binding £2 in case of oscillatory 
responses. Figure 4B shows the clear segregation of the 
parameter sets that give rise to "transient" (green), dam- 
pened (blue) or sustained (purple) oscillatory responses 
in a plane spanned by the receptor-ligand binding rate 
£2 and the speed of I-Smad turn-over (£16 x £17 x 
£19). Oscillations are observed only when £2 is small 
such that ligand binds slowly to its receptor and the 
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Figure 3 Box plots of the parameters corresponding to the different responses types (A) Box plots of parameter sets leading to a 
transient (green) or sustained (red) response. Parameters that differ are mainly k8, /c10, /c1 1, /c13 (shuttling rate from cytoplasm to nucleus, 
formation/dissolution of the Smad dimers, and dephosphorylation of R-Smad), and k5, k6, k14, /c15, kU, k18, k19 (all related to the strengh of the 
feedback). (B) Box plots of parameter sets leading to a transient (green) or oscillatory (blue) response. k16, k\l, k!9 (dynamics of the l-Smad 
mRNA and l-Smad protein) and k2 (binding of TGF-/3 to its receptor) are key determinants of the response kind. Ranges of the uniform sampling 
distributions, as stated in Table S1, are indicated by grey boxes. 



pool of free receptor is depleted gradually (Additional 
file 3, Fig. S6A). As a consequence free receptor is still 
available when I-Smad has downregulated the response 
and ligand can still trigger a further response. Conver- 
sely, if I<2 is large, receptors are rapidly bound to the 
ligand (whose concentration is constant in this model), 
and once the response has terminated, there is no free 
receptor available to induce a new response (Additional 
file 3, Fig. S6B). /c2, and thus the speed with which the 
free receptor concentration decreases, critically deter- 
mine the dampening of oscillations. A 10-fold change in 
the value of I<2 can transform sustained oscillations in 
highly dampened ones (compare panels A and B in 



Additional file 3, Fig. S6). Rapid degradation of I-Smads 
is important for sustained oscillations because otherwise 
all receptors become rapidly sequestered by I-Smad and 
the response is terminated. Accordingly inclusion of 
receptor endocytosis and recycling to the membrane 
(cycling) combined with the removal of the I-Smad 
would allow further oscillatory cycles. 

Each parameter in our simple model integrates the 
effects of many further interactions as may also arise 
from cross-talk. Thus it has been shown that the phos- 
phorylation of Smad in its linker region by Ras-activated 
MAPK induces a cytoplasmic retention of R-Smads [15], 
which in our system would be represented by a lower 
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k8*k11*k13/k10 - Smad deactivation 



k17*k16*k19 - dynamics of l-Smad feedback 
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Figure 4 Impact of kinetic parameters on the type of the TGF-/3 response. (A, C) A strong negative feedback, fast nuclear shuttling of 
Smads and a rapid dissociation of the dinners favour a transient (green) over a sustained response (red). (B, D) Fast production and degradation 
of the l-Smad mRNA and l-Smad protein is required for oscillations to appear, and a low TGF-/3 receptor on-rate enhances oscillatory (blue and 
magenta in the scatter plot, and black to blue in the contour plot) relative to a transient (green in the scatter plot, and green to yellow in the 
contour plot) response. 



shuttling rate into the nucleus (/<8). Interestingly /c8 
indeed strongly influenced the response type. Another 
parameter that appears to be important in determining 
the response characteristics is the binding rate of TGF-/3 
to the receptor (I<2). Our description of the processes at 
the cell membrane is very simple and thus I<2 has also 
to take in account the regulation of TGF-/3 outside the 
cell by soluble sequestering factors and membrane- 
bound co-receptors as well as processes that affect the 
receptor density on the cell membrane. Those auxiliary 
factors play therefore a crucial role in the TGF-/3 path- 
way flexibility. We should stress that all parts of the 
parameter space should be readily reachable for the cell 
and small adjustments in the parameter values should 
thus be sufficient to alter the response type. 

The regulatory impact of cellular protein concentrations 

The kinetic rate constants of a reaction depend on the 
particular protein chemistry. While rate constants may 
be different between species, rate constants are unlikely 
to differ between individuals of one species and even 
more unlikely to differ within a single individual. How- 
ever, during the development of an organism the same 
signaling network can elicit qualitatively different 
responses at different times and locations. We therefore 
wondered whether changes in the protein concentra- 
tions (which can be easily adjusted by an organism or 



result from crosstalks with other signaling pathways) 
would enable the required regulatory flexibility. To find 
parameter ranges that would permit such flexibility we 
repeated our previous screen with different concentra- 
tions of receptors, R-Smad or Co-Smad: for each of the 
three species we first carried out 3 screens where con- 
centrations were increased or decreased from their 
reference concentration c 0 to c 0 /100 or c 0 x 100. We 
then looked for parameter sets that would permit a 
switch between a transient and a sustained output 
response as the protein concentrations changed. Our 
parameter sampling space is huge (23 parameters, with 
most of them sampled over 4 orders of magnitude, 
Additional file 2, Table SI) and a switch could be 
observed for less than 1% of the sets (Figure 5A, black 
bars). When we plotted the parameter ranges for which 
we observed switching we noted that a number of para- 
meter ranges were restricted compared to the initial 
sampling range (Additional file 3, Fig. S7). We therefore 
wondered whether there would be particular parameter 
ranges for which concentration-dependent switching 
would be more frequent. Indeed when we reduced the 
sampling ranges of the parameter values (Additional file 
3, Fig. S7) about 20% of the parameter sets enabled 
switching as the R-Smad concentration was varied, 25% 
as the receptor concentration was varied, and almost 
30% as the Co-Smad concentration was varied (Figure 
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Figure 5 A change in protein concentrations can switch the type of the TGF-/3 response. (A) Percentage of parameter sets that permit a 
switch in the qualitative response (transient versus stustained response) to ligand when TGF-/3 Receptor, R-Smad or Co-Smad concentrations are 
increased or decreased by 100-fold. The parameter sets were drawn from the ranges as specified in Additional file 2, Table S1 (black), or in more 
(grey) and more (light grey) restricted ranges as shown in Additional file 3, Fig. S7. (B) The minimal relative change that is required in the 
concentrations of TGF-/3 Receptor (black), R-Smad (grey), or Co-Smad (light grey) to switch between transient and sustained responses when 
parameters were drawn from the most restricted range (corresponding to the light grey column in panel A and the lightest colour in Additional 
file 3, Fig. S7. 



5 A, grey bars). We notice that the only rates that were 
not restricted while enhancing the fraction of parameter 
sets that permit switching were the rate of ligand- 
TGF/5R unbinding (kl) and of I-Smad mRNA export 
(kl6). These rates thus appear to have very little influ- 
ence on the overall kinetics within the screened range. 
We next wondered what would be the minimal change 
needed in protein concentration to allow the switch. To 
that end we carried out 9 supplementary screens for 
each of the three species where concentrations were 
increased or decreased from their reference concentra- 
tion c 0 over a 100-fold range in multiples of 3, i.e. c 0 (n) 
= c 0 3 n with n = [-4, -3, 3, 4]. Interestingly, while a 
change in the response type was observed most fre- 
quently in response to changes in the Co-Smad concen- 
tration (Figure 5A, grey bars), switches could be 
achieved with much smaller concentration changes 
when the receptor or R-Smad concentration were varied 
(Figure 5B). Thus only a 3-10 fold change in the recep- 
tor and R-Smad concentration was typically required 
while the Co-Smad concentration typically needed to be 
changed by 20-100-fold. The I-Smad Smad6 has indeed 
been reported to inhibit TGF-/3 signaling by sequester- 
ing the Co-Smad Smad4 in an inactive complex [33]. It 
has further been argued that cross-talk between different 
TGF-/3 pathways may be integrated via a competition for 
Co-Smads. Based on our observations such competition 
would need to greatly alter the concentration of avail- 
able Co-Smad to be effective and the receptor and the 
R-Smad would offer a more sensitive point of control. 
Previous models have focused on the dynamical control 
of the TGF-/3 receptor and have shown that this indeed 
offers great regulatory flexibility [29]. Experiments 
further show that the I-Smad may also affect the turn- 



over rate of R-Smads and thus affect their cellular con- 
centration [34]. 

TGF-/J dose-dependent response 

Finally we wondered how different ligand concentra- 
tions would affect the cellular response. The impact of 
different TGF-/3 concentrations have already been stu- 
died by Clarke et al. [35] and Zi et al. [25], but there 
the results were strongly affected by ligand depletion 
since the TGF-/3 concentrations were allowed to go 
down over time because of internalization and degra- 
dation. We were interested how different, but constant, 
stimuli would affect the response - the effect of ligand 
depletion can then be deduced as response to decreas- 
ing ligand concentrations. As we varied the ligand con- 
centration between 0.2 pM and 20 nM we noticed that 
only for a very small fraction of parameter sets (less 
than 0.5%) the response type changed qualitatively as 
the concentration varied. The parameter sets were not 
clustered and a further increase by restricting para- 
meter ranges (as in case of the cellular protein concen- 
trations) could not be achieved. Even though changes 
in the TGF-/3 concentration cannot switch the 
response type in our simulations, the duration of the 
response increases with increasing TGF-/3 concentra- 
tions as previously observed by Zi et al. [25]. This 
increase was, however, insufficient to alter the 
response type according to our definitions. 

The ligand concentration clearly affects the maximal 
response in our simulations, and the transcription factor 
activity increases with the ligand concentration until a 
plateau is reached (Figure 6A). In case of sustained 
responses (but not for transient responses) the peak 
value is reached more quickly at higher ligand 
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Figure 6 TGF-/3 dose-dependent response. (A) The pathway response increases with increasing ligand concentration until a plateau is 
reached. The saturation curve can be described with an exponential function, O peak = max(O pe ^)(1 - exp(-x/r])) where x refers to the ligand 
concentration (which was sampled at 12 concentrations between 0.2 pM and 20 nM), and the parameter 77 is characteristic for the saturation 
concentration. (B, C) A histogram of the distribution of 77, the parameter characteristic for the saturation concentration for (B) the transient 
response, and (C) the sustained response. (D) For the transient set, the parameter that determines if the saturation concentration is high (green) 
or low (dark green) is k2, the binding rate of TGF-/3 to its receptor. (E) For the sustained set, two parameters are crucial. k2 (binding of TGF-/3) 
and k8 (shuttling into the nucleus) set the saturation concentration either to high (red) or low (dark red). 



concentrations (data not shown). The saturation curve 
in Figure 6A can be fitted with an exponential curve, i.e. 
Op e a/< - m & x (0 P eak)(l ' exp (-#//])) where x refers to the 
TGF-/3 ligand concentration (x = 0.2, 2, 6, 10, 15, 20, 50, 
100, 150, 200, 1000, 2000 pM) and the parameter 77 
indicates the concentration range for which the response 
saturates. Histograms of 77 (Figure 6B, C) show that the 
sustained response (Figure 6C) tends to saturate at 
lower TGF-/3 concentrations than transient responses 
(Figure 6B). Moreover, in case of sustained responses 
there is a biphasic distribution in the saturation concen- 
trations with one peak around 0.1 pM and the other 
one around 10 pM (Figure 6C). However, in both transi- 
ent and sustained cases, the transcription factor is able 
to reach similar maximal values (Additional file 3, Fig. 
S8). On the contrary, the maximal output value reached 
by oscillating responses is much lower than in the sus- 
tained and transient case. Our results are mostly in 
agreement with the conclusions drawn by Chung et al. 
(2009), who showed also that transient TGF-/3 responses 
saturate. However, deviating from our results, Chung 
and co-workers observed that also in transient responses 
the peak value is reached more rapidly as the stimulus 
concentration increases. 



For parameter sets that give rise to oscillatory 
responses, changing the input strength and shape does 
not influence the period of oscillation but modulates the 
evolution of the oscillations amplitudes (data not 
shown). When exposed to sustained, high TGF-/3 con- 
centrations the amplitude of oscillations starts to decay 
from the beginning. When the TGF-/3 concentration 
raises progressively, the amplitude of oscillation first 
raises and then decays, reflecting two competing phe- 
nomena: the amplitude of oscillations tends to be pro- 
portional to the input, but at the same time the 
sequestration of the receptor by the inhibitor leads to a 
dampening of the amplitude. 

We next investigated in how far the kinetic para- 
meters can influence the saturation concentration (Fig- 
ure 6B, C) and the maximal output value at saturation 
(Additional file 3, Fig. S8 and Additional file 3, Fig. S9). 
For transient responses it is mainly the rate of ligand- 
receptor binding, /c2, that determines the saturation con- 
centration (Figure 6D and Additional file 3, Fig. S10). In 
case of slow binding higher concentrations of ligand are 
required to saturate the receptors. The saturation con- 
centration for sustained responses are determined both 
by the receptor-ligand binding rate, /c2, and by the 
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cytoplasm-nucleus shuttling rate, /c8 (Figure 6E and 
Additional file 3, Fig. S10). Fast shuttling enables more 
rapid deactivation of Smads as based on observations by 
Hill and coworkers [3] dephosphorylation is restricted 
to the nucleus in our model. As discussed above I<2 and 
kS have both been reported to be modulated by other 
processes. The saturation concentration can therefore 
also be adjusted by cross-talk. 

The different saturation concentrations are likely 
important for the TGF-/3 response as different genes can 
be activated or repressed depending on the nuclear 
Smad complex concentration. While the mechanism by 
which different concentrations of the nuclear transcrip- 
tion factor complex translate into different transcrip- 
tional responses has not been resolved, likely 
mechanisms include promotor selection based on differ- 
ences in the promoter binding-site affinities, cross- 
repression, and the establishment of a reciprocal of 
repressor gene expression [36,37]. 

Proportional "faithful" responses 

When ligands of the TGF-/3 family act as a morphogen, 
as it is for example the case for Dpp in Drosophila or 
Activin in Xenopus, cells must finely sense extracellular 
concentrations and transduce this signal inside the cell. 
We therefore looked for parameter sets leading to a 
response proportional to the input which we term 
"faithful". The parameter sets that gave rise to anything 
but sustained responses (i.e. transient, oscillatory, non- 
responsive, undefined responses) to sustained ligand 
exposure can already be discarded. Those parameter 
sets that gave rise to sustained responses to sustained 
ligand exposure we sought to analyse further with 
dynamic input signals. Here we used as input a function 
that first linearly increased from 0 to 720 pM for 5 



hours and then linearly decreased to zero over the next 
5 hours (Figure 7A). To screen our simulations for 
"faithful" parameter sets we normalized both the input 
and the output with respect to their respective highest 
values, and calculated the squared residuals R between 
input and output according to 

R = E ( in P ut ; " output;) 2 xhe 10% sets with the low- 

j 

est residual were classified as "faithful" and the 10% sets 
with the highest residual were classified as "unfaithful" 
for further analysis (Additional file 3, Fig. Sll). 

A response is faithful if the output is proportional to 
the input over time, i.e. y ou t P ut{t) = a x y input (t), where a 
is the proportionality coefficient. This requires (i) that 
the output adapts rapidly to changes in the input, and 
(ii) that the response does not saturate, i.e. max(y output 
(£)) <max{O pea1 ^ ) which is the case if the proportionality 
coefficient a is low and/or the maximal response value 
max(O peak ) is high. Those requirements are reflected in 
the constraints on the parameter values (Figure 7B-C) 
for faithful responses, i.e. a low binding rate of TGF-/3 
to its receptor and a low phosphorylation rate prevent 
early saturation of the output, while a relative weak 
feedback and a low binding rate of the I-Smad to the 
receptor prevent a premature termination of the 
response. We have previously discussed the regulation 
of the binding rate of TGF-/3 to its receptor, I<2 and 
thus now focus on the feedback. The different I-Smads 
have been shown to vary in their effects. Thus Dad, the 
Drosophila I-Smad, appears to interfere mainly with the 
BMP-like pathways (Tkv and Sax receptor dependent 
pathways) but not the Activin-like Babo-dependent 
pathway [38]. Inhibition by vertebrate Smad6 and 
Smad7 can be achieved by sequestration, enhanced 
degradation, or an impact on phosphorylation. The 




„. ? % 3 10" 4 10" 2 10° 10 2 
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Figure 7 Parameter Dependency of Faithful Responses. (A) To investigate the faithfulness of the response (red curve), we re-analysed those 
parameter sets that had produced sustained responses to a sustained input (red parameter sets in Additional file 3, Fig. S5) with time-varying 
inputs (linearly increasing and then decreasing TGF-/3 input concentration, black curve). Based on the squared residuals (grey area) between the 
normalized inputs (black line) and outputs (red line) we defined faithful and unfaithful responses as those in the first and last 10-quantile 
respectively. (B) An inefficient activation of the TGF-/3 receptor and a weak negative feedback favours faithful (light red to yellow) over unfaithful 
(black to salmon-pink) responses. 
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different processes likely have different efficiencies and 
this will determine the efficiency of the negative 
feedback. 

Our results indicate that under certain parameter 
restrictions the extracellular concentration is directly 
reflected in the output concentration. In that case, TGF- 
P can act as a morphogen, conveying positional infor- 
mation and determining cell-fate, subjected to the set of 
activated and repressed genes. 

Conclusions 

The duration of the signaling response is thought to be 
an important factor influencing the cell's phenotypic 
response to TGF-/3. We have employed a very simple 
model of the TGF-/3 network to better understand the 
mechanistic basis of the observed signaling plasticity. 
We find that the qualitative response (transient, sus- 
tained, oscillations, proportional responses) to a con- 
stant ligand exposure can indeed be changed by altering 
the value of a single parameter value. Since we consider 
a simple model each parameter value represents a wider 
range of processes and our observation thus implies that 
both changes in protein concentration as well as cross- 
talk between signaling pathways can alter the qualitative 
response to a TGF-/3 stimulus. Many more complicated 
models for TGF-/3 signaling as well as for other signal- 
ing networks have been proposed already. To better 
understand the regulatory impact of cross-talk it will be 
important to connect experimentally validated models 
for the TGF-/3 network also to those for other pathway 
models. While many kinetic parameters have been mea- 
sured an important parameter that remains often 
unmeasured is the protein concentrations. To better 
predict the responses in different cell types it will be 
important to obtain quantitative information on protein 
abundance in different cell types - and eventually in 
individual cells. 
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